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Abstract 

A detailed investigation of the low-energy chiral expansion is presented within 
a model truncation of QCD. The truncation allows for a phenomenological 
description of the quark-quark interaction in a framework which maintains 
the global symmetries of QCD and permits a 'i/Nc expansion. The model 
dependence of the chiral coefficients is tested for several forms of the quark- 
quark interaction by varying the form of the running coupling, a{q^), in the 
infrared region. The pattern in the coefficients that arises at tree level is 
consistent with large Nc QCD, and is related to the model truncation. 



I. INTRODUCTION AND SUMMARY 

Phenomenological approaches to quantum cliromodynamics(QCD) continue to provide 
useful intuition into the nature of the strong interaction, and compliment the more direct 
evaluation via lattice techniques. The utility of these treatments is perhaps most apparent 
in the study of chiral observables where lattice calculations are subject to large uncertainties 
due to the extrapolation to light quarks. At low energies, this aspect of QCD is characterized 
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by chiral perturbation theory (^PT) [0-0]. The coefficients of the chiral expansion are input 
parameters to xPT, and their values are determined from experimental observables. These 
coefficients therefore provide a convenient representation of a large body of data relevant to 
low-energy QCD. 

Quantum chromodynamics is formulated in terms of unobserved degrees of freedom - 
quarks and gluons. The presence of these fundamental constituents of hadrons is inferred 
through the analysis of deep inelastic lepton scattering. Nevertheless, the successful appli- 
cation of effective theories such as xPT to a broad range of low-energy strong-interaction 
phenomena suggests that quarks and gluons may be replaced by local effective hadronic 
degrees of freedom in the low-energy domain. This success is largely due to the separation 
in the hadron spectrum between the Goldstone modes and higher mass states. At inter- 
mediate energies it is not clear that such a description remains effective [Q, nor is it clear 
that explicit quark and gluon degrees of freedom are essential. An ideal perspective on this 
problem would be provided if composite hadron fields and their interactions could be mod- 
eled in a manageable form in terms of the point fields of QCD. Functional integral calculus 
formulates this problem as an exercise in changing the variables of integration from quark 
and gluon fields to hadron fields 0. An obvious advantage of this approach is that the 
effective hadron-field interactions retain knowledge of their subhadronic origin. 

The notion that such a change of variables exists for QCD in the low energy domain is 
implicit in the success of the above mentioned hadronic formulations. The explicit oper- 
ation of changing variables allows the underlying dynamics of the microscopic description 
to influence interactions at the macroscopic level. The goal of this "matching" program is 
then to perform the appropriate change of integration variables in the functional integral 
formulation of QCD; 





Significant progress toward this goal has recently been achieved 
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As the local integration variables in ( |1.1|) are identified with the bare hadron fields, their 
effective interactions are simultaneously defined. This process is the result of an expansion 
about the chiral symmetry breaking ground state [^], and an allocation of internal and 
center-of-mass dynamics. The latter is prescribed by the normal-mode expansion of the 
free kinetic operator of the composite particle in a manner analogous to the interaction 
picture of standard quantum field theory The tree-level effective interactions thereby 

obtained occur through a dynamically regulated "constituent-quark" loop and thus reflect 
the underlying description. The low-momentum (gradient) expansion of these tree-level 
nonlocalities produce finite coefficients, and for the Nambu-Goldstone modes is structurally 
consistent with xPT. 

The direct derivation of the chiral coefficients from QCD is presently inaccessible. How- 
ever, they can be derived from a class of chiral invariant quark-based field-theory models of 
QCD which are distinguished by the form of the quark-quark interaction. In this investi- 
gation the sensitivity of the chiral coefficients to the underlying quark-quark interaction is 
tested for a variety of forms to determine their utility in constraining these models. Previous 
work [^|TU| has demonstrated how these techniques can be used to extract the second-order 



coefficients and those at fourth order associated with vr — vr scattering. In the following 
we extend the previous work by calculating the chiral coefficients Li-Lg, and further by 
investigating the sensitivity of these coefficients to the infrared form of the quark-quark 
interaction. The coefficients Lg and Liq are left for a future investigation, however work in 



that direction has been initiated 11 



We find that in general, for momentum-space interactions of the form as{(f)/<f) the 
reproduction of the accepted values of the chiral coefficients requires the running coupling, 
as{q^)i to have a sufficiently large integrated strength to produce dynamical chiral symme- 
try breaking(DxSB), but is otherwise not acutely dependent on the detailed form. This 
constraint is implemented here by fixing the value of the pion decay constant, which de- 
termines the overall scale. The low-momentum strength of as{q^) is then implied by the scale 
at which the infrared phenomenology is matched to the known ultraviolet form. This scale 



is allowed to vary to investigate the sensitivity of the coefficients, and several two-parameter 
models for as{q^) are employed. 

We also find in particular that the coefficients L5 and Lg are most sensitive to the form of 
the interaction in the infrared. These coefficients are primarily responsible for distinguishing 
the pion, kaon, and eta decay constants and providing higher order corrections to their 
masses [Pl,|l2|. The sensitivity of these mass dependent coefficients is an indication that the 
hadron spectrum is playing a role in the determination of the quark-quark interaction. This 
result is consistent with a previous investigation which shows that the convergence radius 
of the chiral expansion in the current quark mass alone strongly depends on the form of the 



quark-quark interaction in the infrared region [|T3|,|T4[. In particular it was found there that 
the running coupling has to be strong in the infrared region in order to obtain convergence 
of the chiral series in the strange quark sector. 

Finally, we find that a pattern in the coefficients emerges at tree level which is consistent 
with the large N^. expansion in QCD and can be traced to approximations that are made 
to QCD here in deriving the low-energy expansion^. In this way the consequences of the 
model assumptions can be directly observed. The present investigation further provides a 
significant reduction in the number of parameters needed to represent the low-energy QCD 
data mentioned above, and thereby deepens our understanding of low-energy QCD. 

The paper is organized as follows. In section II the path from QCD to xPT is explored. 
In Section III the consequences of the approximations made in Section II, along with the 
model dependence of the results, are investigated. Finally conclusions are offered in Section 
IV. 

II. QCD, THE EFFECTIVE QUARK-QUARK INTERACTION, AND xPT 



^The role of the r/o and the associated anomaly are neglected here in considering the A'^^ dependence 
of the chiral coefficients. This question was first addressed in Ref. [|l5| and is reviewed in Ref. p^. 



A. From QCD to the effective quark-quark interaction 



A global color symmetry model(GCM) |^,^ that is based upon an effective quark-quark 
interaction can be defined through a truncation of QCD as follows. The generating functional 
for QCD in the Euclidean metric is given by 



Z[ip, r],r]] = J DqDqDA exp (^—S[q, q, A^] — qipq + r]q + qr] 



(2.1) 



and can be rewritten as 



Z[ilj,rj,ri]=J DqDq exp —J q{^ + ilj)q + rjq + qr] exp^I^ 



(2.2) 



with W[J] given by exp (W^[J]) = / DA exp {-\F^^F^j^ + J^A^). Here f/, rj and il) are 
external source fields. 

The quantity W[J] has an expansion in gluon n-point functions starting at second order; 



W[J] 



Dll{x,y)rAx)Jl{y) + WR[Jl 



(2.3) 



where PVij[J] involves gluon n{> 3)-point functions. It is worth noting that the ra-point 
functions have mass dimension [mass]". One might therefore hope that for low-energy 
hadron physics the low-dimension functions would provide a good description. 

By replacing the quark field variables in by their source derivatives, the generating 

functional of QCD can be written as 



Z[i),r],rj\ = exp Wr 
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^GCM 



(2.4) 



where ZgcmVI^, V,v]= I DqDq exp (-^gca/ [V^, 9, q]+m + QV) with 



SGCM[ip,q,q] = I d xd y<^q{x)[{^^ + ilj{x))6{x-y)\q{y) 

+ ^-lj-{x)D{x~y)jM 



(2.5) 



Here j^{x) = q{x)—^yq{x) is the quark color current, and for convenience a gauge for the 
gluon propagator D^^^{x — y) = 5ah^^j.uD{x — y) is employed. From here forward we work 
within the model truncation defined in 



The primary benefit of tliis truncation is tliat a reasonably solvable model is obtained, 
which is nevertheless sufficiently general to address a variety of phenomenological issues such 
as the role of quark-quark interactions in effective hadronic field theories. This model as 
well maintains the global symmetries of QCD and permits a,l/Nc expansion. 

The primary loss of working at this level is that of the local color gauge invariance of 
QCD. The consequences of this loss are unclear, but are determined by the operation of Wr 
in (|2.4|) . The approximation of a local symmetry by a global symmetry is similar to the 
approximation of general relativity by special relativity. If the relevant field is sufficiently 
weak in the region of interest, then such an approximation is reasonable. In the case of 
localized color-singlet states one might hope that color neutrality could provide such a 



scenario 
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It should be noted that the Nambu-Jona-Lasinio(NJL) model [|T^ is obtained from ( p.5|) 
in the limit D{x — y) = S{x — y)/7Vl^, with Ai the appropriate mass scale. The chiral 



coefficients in the NJL model have been investigated |T8|-|20|], and are a limiting case of 
the present investigation. Our interest here is with the more general question of the model 
dependence of these coefficients. The present description also allows the discussion of higher 
mass excitations due to the nonlocal interaction. 



B. From the quark-quark interaction to xPT 

1. Bosonization and saddle-point expansion 



The meson sector of the variable change implied in ( p, . 1| ) is revealed by first identifying 
field combinations (currents) with the transformation properties of mesons. This is achieved 
through a Fierz reordering of the current-current term of the action ( p.5|) to obtain 
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^f^(x)D{x-y)r,{y) = -l-j\x,y)D{x-y)j\y,x), (2.6) 

where J^{x,y) = q{x)A^q{y) and the minus sign in ( |2.6| ) arises from the Grassmann nature 
of the quark field variables. Here the quantity is the direct product of Dirac, flavor SU{3) 



and color matrices; 

A«= i » (^Ir.^Aj) <8 (ilc.^Aj) . (2.7) 

which contains, in particular, color singlet qq combinations. It should be noted, however, 
that there are also color octet qq combinations present in ( |2.7| ). An alternate color Fierz 
reordering, 

8 4 2 ^ 

X! {K)kl = l^^ilhj + o XI ^rnikemlj, (2.8) 

a=l m=l 

eliminates the color octet qq sector in favor of color triplet-antitriplet qq combinations and 
leads naturally to baryons ||^. This alternate approach, although natural for the investiga- 
tion of baryons, is unnecessary for the investigation of meson interactions of interest here. 
The interested reader is encouraged to consult Ref. and references therein for details of 
the baryon sector. 

Having identified field combinations with the transformation properties of mesons, the 
current-current term of the action ( |2.5| ) is eliminated by multiplying the partition function 
by unity in the Gaussian form 

B%x,y)B%y,xy 



1 = AT J DB exp - J d^xd^ 



2g'D{x-y) J ^^'^^ 
and shifting the bilocal-field integration variables as B^{x, y) B^{x, y) +g'^D{x~y)J^{y, x) 



2T[|. This requires in particular that the bilocal fields B {x,y) display the same symmetry 

transformations as the bilocal currents J^{y,x) |1TT| . 

The partition function now has the form Z[ij] = J\f j DBDqDq e-^I^'^'^-^?] where 

B\x,y)B\y,x) 
2gWix-y) 

(2.10) 



S[^P,q,q,B]= d'xd^yqix) {^^ + ^ix))5ix-y) + A'B%x,y) g(y) + 



The action ( |2.10| ) is quadratic in the quark fields which allows the Grassmann integration 
to be performed by standard methods. The resulting expression for the partition function 
in terms of the bilocal-field integration is Z [■?/']= A/" / DB e^"^'^''^' where the action is given 
by 



S|^,H]^-m4G-]+/.SA^Mf^, (2.11) 

and the quark inverse Green's function, G~^, is defined as 

G-\x, y)^i^. + ^lj{x))S{x-y) + A'B\x, y). (2.12) 

This replacement of the quark-field integration with the bilocal-field integration repre- 
sents an exact functional change of variables. Observable quantities extracted from the 
partition function are unaffected by the variable change, but are now expressed in terms of 
effective (meson) degrees of freedom. A benefit of this is that the effective mesonic interac- 
tions, which are generated from the quark-field determinant in ( p.ll|) , represent a summation 



of quark processes, and are easily exposed by expanding in powers of the bilocal fields. The 
structure of these interactions is illustrated in Fig.l. At this level the bilocal fields interact 
through a bare quark loop as in Fig. la, and do not readily display the dynamics expected 
of quark bound states of QCD. However, as the notion of bare mesons is developed, this 
picture of their interactions is simultaneously refined. 

In anticipation of dynamical chiral symmetry breaking, bare mesons are defined in terms 
of the fluctuations about the saddle point of the action (which is equivalent to the classical 
vacuum). This choice of an expansion point harbors profound dynamical consequences 
in that it largely determines both the structure and interactions of the bare mesons. In 
particular, this choice leads to the rainbow Dyson- Schwinger equation of the quark self 
energy, and the ladder Bethe-Salpeter equation for the internal structure of the bare mesons. 
More importantly, as a result of grouping this particular class of diagrams into bare mesons, 
the expansion about the classical vacuum leads to results for the chiral coefficients which 
are consistent with large Nc QCD, as is discussed in Section III. 

The saddle-point of the action is defined as 4^ =0 and is given by 



Bo, ■4^=0 



Bl{x -y)= g'Dix - y)tT A'Goix - y) . (2.13) 



These configurations are related to nonlocal vacuum condensates and provide self-energy 
dressing of the quarks through the definition S(p) = A^Bq{p) = i ^ [^(p^) — 1]+-B(p^), where 
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and 

This dressing comprises the notion of "constituent" quarks by providing a mass 

Their role as constituents is further displayed by expanding the bilocal fields about the 
saddle point, 

B\x,y) = B',{x-y) + B%x,y), (2.17) 

then examining the effective interactions of the fluctuations, B. These interactions are 
produced by the quark-field determinant TrLn (^^ + S + A^B^^, as is illustrated in Fig. lb. 
There it is seen that the fluctuation-field interactions now occur through the constituent- 
quark loops. 

The connection between the bilocal fluctuation fields and the local fields of standard 
hadronic field-theory phenomenology remains to be shown. The bilocal fields contain in- 
formation about internal excitations of the qq pair in addition to their net collective or 
center-of-mass motion which is to be associated with the usual local field variables. A sep- 
aration of the internal and center-of-mass dynamics is achieved by considering the normal 
modes of the free kinetic operator of the bilocal fields in a manner which is analogous to 
the interaction representation of standard quantum field theory. Details of the localization 



procedure can be found in Refs. [§] and p3|. The process amounts to a projection of the 
bilocal field B^ onto a complete set of internal excitations with the remaining center-of- 
mass degree of freedom represented by the coefficients 0^(-P) = / d'^qB^ {P, q)T^{P, q) . The 
bilocal fluctuations can thus be written as 

B'{P,q)=j:<P:{PK{P,q). (2.18) 



The functions are in general eigenfunctions of the the free kinetic operator of the 
bilocal fields. At the mass shell point, = — M^, they satisfy the homogeneous Bethe- 
Salpeter equation in the ladder approximation for the given quantum numbers 6 and mode 
n. This modal expansion is then used to localize the action. 

At tree level the local fields 0^ interact through a dynamically regulated constituent- 
quark loop, as is illustrated, for example, in Fig. 2. These "effective interactions" thus reflect 
the underlying QCD structure. The intrinsic nonlocality plays a dual role in the subsequent 
description of physical phenomena. First, when sufficiently short length scales are probed 
as in the large momentum behavior of hadronic form factors, the nonlocal structure is 
directly observed ||TT|j2^ . Second, independent of external probes, the nonlocality provides 
a regulation of internal loop integrations, and serves to suppress hadron-loop effects . The 
present approach can also accommodate extensions of low-energy effective theories through 
the consideration of the higher mass states, and therefore provides a consistent framework 
in which many of the issues facing hadronic field theories might be addressed. 



2. Dynamical chiral symmetry breaking 

In the following discussion, D^SB is associated with the occurance of a massless Gold- 
stone mode that is related to the dynamical generation of a scalar amplitude in the quark 
self energy in the limit of vanishing quark mass. We begin by considering the axial-vector 



Ward identity in the chiral limit given by |]26| , |27 | 



m=0 



G-\q + P/2)j, + j,G-\q-P/2). 



(2.19) 



It is well known 26,27 that in the chiral limit the axial- vector vertex contains a zero- 



momentum pole of the form 



(2.20) 



associated with the massless Goldstone mode. It should be noted that in ( ^.2(J| ) the quark- 
pseudoscalar vertex F5 is also evaluated in the vicinity of P = 0, and is therefore a solution 



10 



of the homogeneous Bethe-Salpeter equation for the pseudoscalar bound state. 
Operating on ( p.20| ) with and comparing with (|2.19| ) obtains 

Big') 



r5(0,g) = 275 



(2.21) 



This is the Goldberger-Treiman relation for the quark-pseudoscalar vertex. The fact that the 
quark self-energy function B occurs as the residue of the zero momentum pole in the quark- 
axial- vector vertex is equivalent to a statement of Goldstone's theorem in this context. It is 
also readily verified that the ladder Bethe-Salpeter equation for the pseudo-scalar Goldstone 
mode reduces to the self-energy equation (|2.15|) ||26| . |27|] . 

Since our interest here is the effective action for the Goldstone modes, we neglect all of 
the higher mass fluctuations present in the bilocal fields. This implies that the full bilocal 
field of Eq.( |2.17| ) can be written, using the expressions (|2.13|) and ( p.21|) for the saddle point 
and the Bethe-Salpeter amplitude of the Goldstone modes respectively, as 



A^B\x, y) = E{x -y) + B{x - y) 



X + y 



- 1 



(2.22) 



where U^^x) = PfiU {x) + PlW (x) with Pji l the standard right-left projection operators. For 
the SU{3) flavor case under consideration here the chiral field U is defined as U = e^^"'^"^^^ . 
It should be stressed at this point that we have not integrated over the higher mass states, 
but have simply neglected them. The effect of including and integrating over the higher 
mass states is addressed in Section III. 



3. The low-energy expansion 

For the application to low-energy observables, an expansion of the action to fourth order 
is now considered. The usual chiral power counting is observed p|,p!2|. In order to preserve the 
chiral invariance of the full action ( |2.5| ), the quarks have to be coupled to the external source 
field iplx), which transforms in a certain way under chiral rotations |1||13- -'■^ performing 
the gradient expansion it is important to keep the x dependence of this external field. After 
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carrying out the gradient expansion to fourth order, we will employ the equation of motion 
which is obtained at second order and depends on the external field iplx), and then finally we 
will identify iIj{x) with the current quark mass matrix. Failure to keep the x dependence of 
ip to the very end violates chiral invariance and will render unphysical results for some of the 
low-energy coefficients. This approach differs somewhat from the previous work of Refs. 
and [|l^] where the equation of motion is not employed. However, there the mass- dependent 



fourth-order coefficients are not considered. 

We consider here only the real contribution to the effective action. The imaginary con- 
tribution, which contains the Wess-Zumino term, has also been investigated in Refs. and 
10| , and the interested reader is encouraged to consult these references for more details. 
The restriction of the fluctuations to Goldstone modes with UW = 1, as in ( p. 221) , 



entails that the second term of the action in Eq. (|2.11| ) is independent of the flelds U and 
can therefore be neglected. The real contribution to the action is then given by 



1. 



2 

where is, from ( |2.12| ) and ( |2.22D , given by 



S = Re[S] = --TrLn G'^ G'^ , (2.23) 



G-\x, = 7 ■ d.A{x -y) + ^ 6{x - y) + B{x - y)U, . (2.24) 



By expanding the logarithm and dropping irrelevant constant terms, Eq. (|2.23| ) can further 
be written as 

S = - Y^^Tiia + h + c + df , (2.25) 

where a, 6, c, and d are non-commuting operators formed from A, S, f/5, and and are at 
least of order one, one, two, and three in chiral counting respectively. The explicit form of 
these operators is given in the appendix. 

The effective chiral action to the desired order is now obtained by expanding the sum in 
Eq.( p^.25D and expanding the operators a-d in gradients. The result to fourth order is (in 
Euclidean space) 
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s 



d X 



J w 



tr 



J w 



L, (tr [{d,U){d,U^) 



Lot! 



{d,U){dM^)] - tr \{d,U)idM^) 



(2.26) 



{d,U)id,U^)idM){dM^) 



{d,U){d,U^)iUx^ + xU^) 



Lsti IxU^xU^ + Ux^Ux^ 



where x{^) = ~2(gg)'?/'(x)//^ and the remaining trace is over flavor. In obtaining this result 
the equation of motion 



{d'U)U^ + {d,U)id,U^) + ^{xU^ - Ux^) = 



and the SU (3) relation [jl2 



tr 



id,U)idM^)id,U)idM^)] = - (tr [id,U)id,U^) 



-tr 



tr 



2tr 



have been used. Explicit forms of the coefficients are given in the appendix. 



III. RESULTS AND DISCUSSION 



(2.27) 



{d,U){d,U^)idM){duU^)] (2.28) 



Several conclusions can be drawn directly from the low-energy expansion (|2.26| ). It is 
immediately apparent that the coefficients L4, Lg, and Lj vanish. It is also evident, by 
application of the SU{3) relation ( p.28| ) (see appendix), that L2 = 2Li. These relationships 
are expected in the large limitQ of QCD [|T^. The fact that they are produced here is 
perhaps not too surprising and can be linked to our truncation of the QCD action to include 
only the gluon two-point function. 

With only a two-point quark-quark interaction, the large Nf. limit leads to a description 
of mesons as a sum of ladder exchanges. Our description of mesons as fluctuations about 



^In the presence of the Ua{^) anomaly the coefficient L7 is of order A'^^ |15,12|. Our neglect of 
the rjQ here leads to the vanishing of Lj. 
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the saddle point of the action is equivalent to the 1/Nc expansion and, in this model trunca- 
tion, leads directly to the ladder approximation. Our further neglect of higher mass states, 
explicitly excludes intermediate states of pure glue which are 'Wc suppressed" . A departure 
from this tree-level pattern in the coefficients would therefore have to arise in the present 
formalism by including and integrating over the higher mass mesons, which we have explic- 
itly excluded in Eq. (|2.22| ). The role of the underlying description is thus clearly displayed 
in the pattern of the chiral coefficients. 

Examples of the diagrams that are generated by integrating over higher mass mesons 
are illustrated in Fig.3. The diagram of Fig.3a is of order one in Nc counting and produces 
departures from the large N^. relations, while the diagram of Fig. 3b is of order N^. and 
produces, for example, the p-pole in vr-vr scattering. All of the contributions that we are 
presently considering are of order N^. and arise from a single quark loop'^. 

The remaining nonzero coefficients must be evaluated numerically. These depend explic- 
itly of the values of the self-energy functions A and B in Eqs.( |2.1^ ) and ( p.l5| ) respectively, 
and are therefore implicitly dependent on the quark-quark interaction D. The procedure 
is then to select a form for the function solve the coupled nonlinear equations ( [^.14[ ) 
and ( [^.15D for A and B respectively, and then evaluate the pion decay constant /^r, the 
condensate (gg), and the coefficients Li, L3, L5, and Lg. 

The quark-quark interaction D has the form 



where s = g^, and we investigate three different two-parameter models for a(s); 




4A2 ln(s/A2 + e) 



a2{s) = Tcd 



s^ + A ln(s/A2 + e) 



(3.2) 



^The double trace terms proportional to Li and L2 in ( p. 26 ) arise here from the SU{3) relation 
( 2.2§| ) and originate from a single quark loop. 
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a3{s) = nd 



1 + xe~"^^ 
ln(s/A2 + e) 



Each of these forms incorporates the one-loop perturbative result for large s (here A = 
0.2GeV and d = 12/(33 — 2Np) = 12/27), and extrapolates differently into the low- 
momentum region. The two low- momentum parameters, x ^-^^d A, are varied with the 
pion decay constant held fixed at /,r = 86MeV. This value is appropriate at zero-momentum 
rather than the pion-mass-shell value of 93MeV, however the results are not very sensitive 
to this small difference. By fixing the overall scale of D^SB is fixed. The remaining 
independent parameter is associated with the matching scale to the perturbative form. 

The running coupling for the three cases listed in Eg. ( |3.2| ) are plotted in Figs. |, |^, 
and ^ along with the corresponding solutions of Eqs.( p^.l3 ) and ( p^.l5| ) for the self-energy 



functions. In all three cases as the matching point to the perturbative form is decreased to 
lower momentum, the infrared strength must be increased to maintain the fixed value of /^r. 
Thus the integrated strength of a is largely constant. This trend is also present to a lesser 
extent in the self-energy functions. 

The first model, ai in Eq. (|3.2|) , has been used in previous investigations of the present 
type ||2^. There the parameters were fixed at A = 0.002GeV^ and x = 1.14GeV, which 



leads to a slightly lower value of /^r. The small A limit of this model obtains a matching 
point near zero momentum and a delta-function behavior in the quark-quark interaction D. 
This limit has been used previously to model confinement |^^. The infrared contribution 
to the second model, a2, generates a l/g*^ singularity in the quark-quark interaction D in 
the limit as A — 0. Such a singularity has also been considered previously as a model of 
confinement ||30|. The 1/g^ form falls much slower than the Gaussian in ai, and hence leads 



to much higher matching scales. Finally the third model as has been chosen here to be 
structurally different from ai and 0:2 in order to further illustrate the independence of the 
results to the details of the low-momentum parameterization. The corresponding results for 
the low energy coefficients Li, L3, L5, and Lg are displayed in Tables I-III, respectively. 
In all of the three cases the same pattern is observed: The coefficients Li and L3, which 
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are responsible for vr-vr and K-K scattering, are nearly independent of the form of a{s) and 
therefore on the form of the quark-quark interaction, provided that the integrated strength of 
a{s) is fixed by /^r. On the other hand, the mass dependent coefficients, L5 and Lg, are more 
strongly dependent on the actual form of the interaction. For example with L5, in order to 
reproduce the experimental value, forms of a{s) with a small matching scale, i.e. which are 
relatively strong in the infrared region, are required. This observation is in coincidence with 
the result of Ref . fl^ , where it is shown that quark-quark interactions with a low matching 



scale are also required to achieve convergence of the chiral series in the strange quark sector. 
Furthermore an explanation for the success of the "delta-function-plus-tail" type models 



IS 



(obtained for example from ai in the limit A — 0) in describing chiral observables p9 
offered by this fact. 

We also find that the results for the fourth-order coefficients are rather insensitive to the 
asymptotic UV tale of even omitting this tail completely gives no significant changes, 
again provided that is fixed. 

An increased accuracy in the experimental determination of the coefficients would make 
tighter restrictions on the quark-quark interaction, however the additional investigation of 
higher mass excitations is clearly required to gain detailed information on its infrared form 
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It has frequently been stated with regard to the fourth order coefficients that QCD 
"seems to predict that deviations from the lowest order chiral relation must be in such a 



form as to reproduce the low energy tails of the light resonances, in particular the p." [12 



Here we have explicitly neglected in Eq.( p.22|) the qq fluctuation associated with the p, and 



illustrated in Fig. 3 how the p-pole contribution would arise. One might then ask: What is 
the mechanism that produces the p-tail-like contribution to the coefficients here? 

This question is easily answered by again considering the diagram of Fig.2. There is a 
qq pair in the intermediate state which arises from the quark loop structure of the interac- 
tions. The integrands of these quark loops are peaked at a momentum qpeak such that the 
constituent mass of Eg. ( p. 161) gives 
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M(g2) ^ M(g2^,,) ^ 300 - 400MeV. (3.3) 

This qq pair can have the quantum numbers of the p, and carries sufficient mass to contribute 
the p-tail effect away from the pole. 

IV. CONCLUSIONS 

We have made a detailed examination of the low-energy chiral expansion from the stand- 
point of the model truncation of QCD given in Eq. (|2.5|) . The structure of the model main- 
tains the global symmetries of QCD (including global color symmetry), and permits a,l/Nc 
expansion. The infrared momentum dependence of the quark-quark interaction is phe- 
nomenological input to the model; here three different two-parameter forms are investigated. 

We find that by truncating QCD to include only a two-point quark-quark interaction 
and describing mesons as fluctuations about the saddle point of the effective action, one 
obtains a pattern in the chiral coefficients which is consistent with large N^. results in QCD. 
This conclusion can be understood by considering the l/N^. expansion within the model 
truncation, and provides a direct link between the model assumptions and consequences for 
physical observables in QCD, independent of the phenomenological treatment of the quark- 
quark interaction. The structure of the underlying theory is in this way displayed by the 
pattern in the chiral coefficients. The departure from the large result is provided here 
by the integration over higher mass states. 

We find that by fixing the pion decay constant, an integrated strength of the running 
coupling is prescribed. This sets the scale for D^SB. The remaining independent parameter 
is associated with the matching scale to the perturbative form of the running coupling. The 
chiral coefficients Li and L3, which are related to tt-tt and K-K scattering data are nearly 
insensitive to this scale. It appears, therefore, that any chiral quark-quark interaction which 
is capable of D^SB can be expected to reproduce these coefficients and the corresponding low 
energy meson scattering data. However, some constraint on the matching scale is provided 
by the sensitivity of the mass-dependent coefficients L5 and Lg, which favor interaction forms 
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that are strong in the infrared domain. 

Finally we conclude that the model truncation that is employed here reproduces low- 
energy QCD, as represented by xPT, quite well. More importantly, this model is not limited 
to low energy and might therefore be used to extend low-energy effective theories. 
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APPENDIX A: 



The operators in Eq( p.25[ ) are defined as 



a 

X ' 

b = IbP^uI - 1]B 



c = -[BU^^^ + tpU^B + ^V] (Al) 



X 



where for example 



< xi\Af,\x2 >= d^^^A{xi - X2), < Pi\A^\p2 >= ipi^A{pl)5{pi - P2) 

< Xi\B\x2 >= B{xi - X2), < Pi\B\p2 >= B{p\)5{pi - P2) (A2) 
< x,\BU,\x2 >= B{x, - X2)U, (^4^) ' < ^il^^5|P2 >= ^5 (^^) U,{p, - P2) 



2 ; ' ^^'^ (27r)2 V 2 

Prom the discussion presented in the text, one can then obtain the coefficients 

f^^Pj dss(^i^^^\A'' + sAA' + s\A'y + s{B'Y] 

BB' + ^[{B'f + BB"]\ 



{m)lGeV ^ ~^ J 



X 
X 















L3 = 


-(A2- 


2A3 + A1) 


^5 = 


A4 — Ae 




^8 = 


-(A5- 


^A, + ^Ae) 



where 



Ai = cis(Aii + A12 + Ai3 + Ai4 + A15) 



with 



All = ( B'^ + BB" + ^sBB'" + sB'B" 

32 V 3 

= -85(i?i?0^ 4 - KBB'f-^ + 2^[55'" - B'B"] 

64 y x"* x^ 

+2^[3BB" - {Bf]^ 

Fs ( 

Ai3 = -— 3[(5')' + BB"\ + 3s[55'" + 3S'5"] + -\BB"" + 4S'S'" + 3(5")'] 



Ai5 = l^sS' 
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ds{X2i + A22 + A23 + A24 + A25 + A26 + A27 + A28 + A29 + A210) 



+sBB' 
F 



BB'f + sBB'[{B'f + BB"] + ^[{B'f + BB'f 



{B'f + BB" + -BB'" + sB'B" 



3 ysy X J 4:^ ' x^ 

Fs^ /BB'Y 
2Fs 



X 



x^ 



B^B'V 



-Fs 



B^B'V (s^ 



-AA' + -{A'f 



-^^(a'^BB' + ^A^B'f + BB"] + sBB'[AA' + s{A'f] 



+^s^[AA' + s{A'f][{B')^ + BB"] 



Fs^ B^B' 

4 x^ 



(A') 



/\2 



sF B ( 4 ^B^B'x' 



/\2 



X 



+ SS"] + + sAA' + A"] 



2 X 
sF /B\* 



^(^) (a^ + 2sA'A' + -s^AA'Y + -s^A(A'f + -s\A')A 
A \x J \ 3 3 3 / 



A3 — ds(Xsi + A32 + A33 + A34 + A35 + A36 + A37 + A38 + A39 + A310) 



(A8) 



{&{jii^r^BB"]^^lBB' 
64 I 3 ^ ' x^ 



{B'Y + BB" + -BB" 



X 



sB'B" 
BB' 



8{BB'f^ + 2[BB"+{B'f 



X'' 



x^ 
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2 ,BB"'-B'B"\ 
2s'F ( 1 (BB'\' s {B'f + BB "' 

+ o(^^) -3 



?> \A \ X ) 8' ' 

2 D/ \ 2 / „2 „3 \ 



B'B'Y (s' . ., s 
x^ J 
sFB^ /s2 



s^F B^B' 



4 x3 V 3 



4 ,t3 

sFBf 2,B^'B'x'.^,., , ,,,, «o2o/.2^' , 2 oor.oA2 , + 



■s 



[s{A'y + A' A] - -B^B'A^=- + -s'BliB'f + SS" 



4 Y 3 2 3 X 

+ '-A\B'f + SS"]- + ^[s\A'f + sAA' + A^\\ 
2 x X / 



sF fB 



( ^\ (-A^ - 2sA'A' - -s^AA'f + -s^A{A'f + -s^A'A 
\x J \ 3 3 3/ 



I ( -)\s{B'f + A^ + sAA! + s'(A)''] 



(A9) 



As = ^ / o?s(A6i + A62 + A63) (All) 



X^,^-tl^^[A'^sAA'-rs\A'f] 
4 

_sFB' + |S" 

^ 8 x 
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and finally 

F.^ and B..-<|. (A12) 

In the above expressions, the arguments of the functions are s — (f and the prime 
indicates differentiation with respect to s. 
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FIGURES 

FIG. 1. The effective interactions obtained from the quark determinant are illustrated both 

before and after the saddle-point expansion. 

FIG. 2. An example of the effective interactions between the localized mesons is shown. The 
quark lines and vertices are dressed in the rainbow and ladder approximations respectively. 

FIG. 3. Examples of the diagrams generated by integrating over higher mass mesons are shown. 
The diagram in (a) is of order one in Nc counting while that of (b) is of order Nc- 

FIG. 4. The running coupling ai, and quark self-energy functions A and B as shown versus 
s = (f. The parameter choices maintain = 86MeV. 

FIG. 5. Same as Fig.4 using 02- 

FIG. 6. Same as Fig. 4 using 03. 
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TABLES 



a,(s) = 3TTSx'c-^'^/(i^-) + TTd/ln(s/A^ + e) 


A (GeV2) 


X (GeV) 


-{qq)^'^ (MeV) 


Li(0.7±0.5) 


L3(-3.6±1.3) 


L5(1.4±0.5) 


L8(0.9±0.3) 


0.002 


1.4 


150 


0.84 


-4.4 


1.0 


0.88 


0.02 


1.5 


160 


0.82 


-4.4 


1.14 


0.84 


0.2 


1.65 


173 


0.81 


-4.0 


1.66 


0.83 


0.4 


1.84 


177 


0.80 


-3.8 


2.0 


0.87 



TABLE I. The chiral coefficients, calculated using ai, are displayed. The parameter choices 
listed maintain /^r = 86MeV. The "experimental" values in parenthesis are taken from Ref . [3] . The 
quark condensate {qq) is evaluated at 1 GeV. 



oi2{s) = nd{sx^/{s^ + A) + l/ln(s/A2 + e)) 


A (GeV^) 


X (GeV) 


-(gg)i/3 (MeV) 


Li(0.7±0.5) 


L3(-3.6±1.3) 


L5(1.4±0.5) 


L8(0.9±0.3) 


10-^ 


0.83 


162 


0.82 


-4.4 


1.28 


0.87 


10-^ 


1.02 


167 


0.81 


-4.2 


1.60 


0.91 


10-1 


1.83 


173 


0.79 


-3.8 


2.36 


1.00 


1 


2.73 


173 


0.79 


-3.5 


3.0 


1.17 



TABLE IL Same as Table I using 0:2. 
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asis) = 7rd(l + xe-^/^)/ln(s/A2 + e) 


A (GeV^) 


Y 
A 


-(Qa)^/3 (MeV) 


Li(0.7±0.5) 


L3(-3.6±1.3) 


L5(1.4±0.5) 


L8(0.9±0.3) 


0.1 


61.0 


163 


0.82 


-4.3 


1.22 


0.83 


0.4 


24.0 


169 


0.81 


-4.2 


1.48 


0.84 


1.0 


15.3 


171 


0.80 


-4.1 


1.73 


0.88 


2.0 


12.2 


172 


0.80 


-1.0 


1.97 


0.95 



TABLE III. Same as Table I using 0:3. 
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(a) Running Coupling Constant a(s)=37rsxV^^/(4A^)+7rd/In(s/A^ +e) 



Q;s(s)=7rd/ln(s/A +e) 
A=0.4(GeV ), x=l-84(GeV) 
A=0.2(GeV\ x=l-65(GeV) 
A=0.02(GeV^), x=1.53(GeV) 
A=0.002(GeV^), x=l-41(GeV) 
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